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Abstract. In this article, we study the 1 + 3 dimensional relativistic Eulcr 
equations on a pre-specified conformally flat expanding spacetime background 
with spatial slices that are diff'eomorphic to K.^ . We assume that the fluid ver- 
ifles the equation of state p = c^p, where < Cs < ^1/3 is the speed of 
sound. We also assume that the inverse of the scale factor associated to the 
expanding spacetime metric verifies a Cs— dependent timc-intcgrability condi- 
tion. Under these assumptions, we use the vectorfleld energy method to prove 
that an explicit family of physically motivated, spatially homogeneous, and 
spatially isotropic fluid solutions is globally future-stable under small pertur- 
bations of their initial conditions. The explicit solutions corresponding to each 
scale factor are analogs of the well-known spatially flat Friedmann-Lemaitre- 
Robertson- Walker family. Our nonlinear analysis, which exploits dissipative 
terms generated by the expansion, shows that the perturbed solutions exist 
for all future times and remain close to the explicit solutions. This work is an 
extension of previous results, which showed that an analogous stability result 
holds when the spacetime is exponentially expanding. In the case of the radi- 
ation equation of state p = (l/3)p, we also show that if the time-integrability 
condition for the inverse of the scale factor fails to hold, then the explicit fluid 
solutions are unstable. More precisely, we show that arbitrarily small smooth 
perturbations of the explicit solutions' initial conditions can launch perturbed 
solutions that form shocks in finite time. The shock formation proof is based 
on the conformal invariance of the relativistic Euler equations when = 1/3, 
which allows for a reduction to a well-known result of Christodoulou. 
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1. Introduction and Summary of Main Results 

In this article, we study the future-global behavior of small perturbations of a 
family of physically motivated explicit solutions to the 1 + 3 dimensional relativistic 
Euler equations. We assume throughout that the fluid equation of state is 

(1-1) P-c^p, 

where p > denotes the pressure, p > denotes the energy density, and the 
constant > is known as the speed of sound. The explicit solutions model a 
spatially homogeneous isotropic fluid of strictly positive energy density evolving 
in a pre-specified spacetime-with-boundary ([l,c») x M^,g) which is undergoing 
expansion. Here, g is an expanding Lorentzian metric of the form (|1.2p . As is dis- 
cussed in Section [1.1 1 in more detail, such solutions play a central role in cosmology, 
where much of the "normal matter'Q content of the universe is often assumed to 
be effectively modeled by a fluid with the aforementioned properties. 

Although all speeds verifying < < 1 are studied in the cosmology literature, 
the cases = and = 1/3 are of special significance. The former is known as 
the case of the "pressureless dust," and is often assumed to be a good model for 
the normal matter in the present-day universe. The latter is known as the case of 
"pure radiation," and is often assumed to be a good model of matter in the early 
universe. For an introductory discussion of these issues, see jWal84| Section 5.2]. 
As we shall see, in order to prove our main results, we will assume that G [0, 1/3]. 
We simultaneously analyze the cases € (0, 1/3) using the same techniques, while 
the endpoint cases = and = 1/3 require some modifications; when = 0, 
p loses one degree of differentiability relative to the other cases, while the analysis 
of the case = 1/3 is based on the fact that the relativistic Euler equations 
are conformally invariant in this case (see Section |9]). We will also provide some 
heuristic evidence of fiuid instability when > 1/3; (see Section ri.2.3p . 

We assume throughout that the pre-specified spacetime metric g is of the form 

3 

(1.2) 5 = -dt2 + e2^^(*)^(dx^)2, 1](1)=0, 

where e^*-*-* > is known as the scale factor, t e [l,oo) is a time coordinate, and 
(cc^, x^, x^) are standard coordinates on R^. A major goal of this article is to study 
the influence of the scale factor on the long-time behavior of the fluid. Our analysis 
will heavily depend on the assumptions we make on e'^'-*-', which are stated later in 
this section. Metrics of the form ()1.2p are of particular importance because exper- 
imental evidence indicates that universe is expanding and approximately spatially 
flat |EBW92j . According to Steven Weinberg [WeiOSl pg. 1], "Almost all of mod- 
ern cosmology is based on [metrics of the form (jl.2p ]." The physical significance of 
these metrics is discussed more fully in Section [Ol 

Under equations of state of the form p = cj.p, the relativistic Euler equations in 
the spacetime-with-boundary ([l,cxj) x M'^,(7) are the four equations 



In cosmological models, "normal matter" is distinguished from "dark matter;" the latter 
interacts with normal matter indirectly, i.e., only through its gravitational influence. 
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(1.3) D^T"^^Q, (/^ = 0,1, 2, 3), 
where 

(1.4) T^'" = {p + p)u^'u'' +p{g-^r'' 

is the energy-momentum tensor of a perfect fluid, D is the Levi-Civita connection 
of g, and u^^ is the four-velocity, a future-directed {u^ > 0) vectorfield normahzed 

by 

(1.5) gc.pu'^u^ = -I. 

Equivalently, the following equation for holds in our coordinate system: 

(1.6) W° = (l+,9a6wV)V2. 

Under the above assumptions, it is well-known (see e.g. |Spell[ ) that when 
p > 0, the relativistic Euler equations ()1.3p can be written in the equivalent form 

(1.7a) u'^Dc,lnp+ {l + cl)DaU°' =0, 

(1.7b) u^'Dc.u'' + —^U"^Do.\np = 0, (/^ = 0, 1, 2, 3), 

l + cj 

(1.7c) n^"^ ='^^^^^'^ + (.g-l)''^ (//, = 0, 1, 2, 3). 

The tensorfield 11'"^ projects onto the g— orthogonal complement of u^. We remark 
that the system (|1.7ap - (|1.7bp is consistently overdetermined in the sense that the 
component of (jl.7bp is a consequence of the remaining (including p.6p ) equations. 

Other important experimentally determined facts are: on large scales, the matter 
content of the universe appears to be approximately spatially homogeneous (see e.g. 
[YBPSOSQ and isotropic (see e.g. |WLR990 . These conditions can be modeled by 
the following explicit homogcneous/isotropic fluid "background" solution: 



(1.8) p"='>e-3(i+c=)a(*)^ u^'^'s^, 

where p > is a constant. It follows easily from the discussion in Section [3T] that 
()1.8|) are solutions to (|1.7ap - (|1.7bp (where {t,x^,x'^,x^) are coordinates such that 
the metric g has the form (|1.2I) ): see Remark 1 3. II 

The main goal of this article is to understand the future-global dynamics of 
solutions to (|1.7al) - (jl.7bl) that are launched by small perturbations of the initial 
data corresponding the explicit solutions ()1.8|) . In order to prove our main stability 
results, we make the following assumptions on the scale factor. 

Assumptions on 6*^^*^ 
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A3 < 



Al e^^') e C^([l,oo), [1,00)) and increases without bound as t — >■ 00 

J-^ g-20(s) <. ^2 ^ 0^ 

A2 { fZi e""(')^(f^(s)) ds < 00, < c2 < 1/3, 
/J^^ e-"(^) d5 < 00, = 1/3 

.^(•) € C^([l,oo), (0,00)) and increases without bound, < < 1/3, 
< (1 - 3c2) for ah large f], < < 1/3, 

. /^^i 4l) < ^' < < 1/3. 

We now roughly summarize our main results. See Theorem 18.11 and Corollary 
for precise statements. 
Main Results. Assume that the scale factor of g verifies assumptions 
lAll - [A3] and that the equation of state p = c^p holds, where < < 
1/3. Then the explicit solution (p,u^^) defined in (|1.8p is globally future- 
stable. More precisely, if the initial data corresponding to the explicit 
solution are perturbed by a sufficiently small element of a suitable Sobolev 
space, then the corresponding solution of (|1.7al) - (jl.7bp exists classically 
for {t,x^,x'^,x^) e [1,00) X R'^ and remains globally close to the explicit 
solution. In particular, no future-shocks form in such solutions. 

In contrast, if = 1/3, assumption lAll holds, but /J^j^ e"^^^-* rfs = 
00, then the solution (p, u^) is nonlinearly unstable. In particular, there 
exists an open family of initial data containing arbitrarily small smooth 
perturbations of the explicit data whose corresponding perturbed solutions 
develop shock singularities in finite time. 

Remark 1.1. Since we prove our stability results through the use of energy meth- 
ods (which are stable), the precise form of the metric (|1.2p is not essential and was 
chosen only out of convenience; our stability results would also hold under suitably 
small perturbations of g (belonging e.g. to a Sobolev-type space of sufficiently high 
order). However, it is not clear whether or not the shock formation result is stable 
under perturbations of the metric; see Remark 19.21 



Remark 1.2. As we will see in Proposition 13. II and Corollary [321 (see also Section 
ll.2.3p . spacetime expansion generates dissipative terms in the relativistic Euler 
equations. This dissipative mechanism is the main reason we are able to prove 
our future stability results. A related fact is that our proof of future stability 
also applies to the spacetime-with-boundary {g,[l,oo) x T'^). We remark that a 
framework for proving small-data global existence for quasilinear hyperbolic PDEs 
in the presence of strong dissipation has been developed in e.g. |Hsi97| (sec also 
[NisTSj for the case of quasilinear wave equations), but its assumptions roughly 
coincide with the dissipative effect generated by e'^'-*^ = e^* (where iJ > is a 
constant) and are therefore much stronger than[AT]-IA3l 

Some important examples of scale factors that appear in the cosmology literature 
and that verify lAll - IA3I include i) the case of exponential expansion, in which 
e^(*) ~ e-^* for some "Hubble" constant H > (choose ^(fJ) = e'^ for some small 
constant q > when < < 1/3). Future stability in this case is implied by the 
results proved in [RS09| . |Spell[ under the assumption < < 1/3. ii) the case 
of accelerated power-law expansion, in which e^'*^ ^ t'^ for some constant Q > 1 
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(choose I^{n) — e''^^ for some small constant g > when < < 1/3). The 
significance of these expansion rates will be discussed in more detail in Section [01 

1.1. Motivation and connections to previous work. We are interested in the 
system p. Tap - (jl.7b[) because it plays a fundamental role in the standard model 
of cosmology. More precisely, physicists/cosmologists often couple the Einstein 
equations of general relativity to the relativistic Eulcr equations and sometimes also 
to additional matter models. They then study explicit solutions of the resulting 
coupled system and make predictions about the long-time behavior of the universe 
based on the properties of the explicit solutions. Now a basic requirement for 
the explicit solutions to have any long-time predictive value is that they he future- 
stable under small perturbations of their data] this basic requirement is the primary 
motivation behind our investigation. 

In the coupled problem mentioned above, the spacetime metric g is not pre- 
specified, but is instead one of the dynamic variables. In this article, we don't 
address the coupled problem, but rather the behavior of the fluid matter under the 
influence of a large class of pre-specified spacetime metrics g. Our hypotheses I Al I - 
IA3l on g are meant to roughly capture the main features of the kinds of metrics that 
may arise in a large class of physically motivated expanding solutions to the coupled 
problem. Readers may consult e.g. |Ren05b] . |Ren06j for background information 
on the coupled Einstein-matter problem in the presence of expansion. The physical 
relevance of these expanding metrics is the following: experimental evidence sug- 
gests that our own universe seems to be undergoing accelerated expansion. There 
are many experimental references available on the issue of accelerated spacetime ex- 
pansion; !RFC+98j , [PR99] are two examples. However, as noted in e.g. |Ren05aj . 
|Ren06| . |Rin09| . the precise rate and mechanism of the expansion of our 
universe have yet to be determined. This uncertainty has motivated us to 
make the very general hypotheses lAll - IA31 We expect that these hypotheses will 
be roughly consistent with the behavior of most of the physically relevant metrics 
g that arise in the study of accelerated expansion in the coupled problem. Now 
for metrics of the form (|1.2p . accelerated expansion (i.e. JjC^^*-', ^e^^*-' > 0) is 
equivalent to (where Ljj{t) ■^Vl{t)) 

(1) w > 

(2) f^uj + uj'^ > 0. 

We note that these conditions are not equivalent to lAll - [A3l it is not difficult to see 
that (1) — (2) may hold even if lA2l - [A3] fail (for all functions ^i-)), and conversely 
lAll - IA3I may hold without the expansion being accelerated for all time (we don't 
even need cj to be differentiable in order to prove our results). Nevertheless, many 
of the scale factors appearing in the cosmology literature verify (1) — (2) and also 
IA1I - IA3] Later in this section we will further discuss the two previously mentioned 
examples of exponential expansion and accelerated power-law expansion. 

Our main results are further motivated by the following: based on our previous 
experience |RS09| . |Spell| , in certain accelerated expanding regimes, we expect that 
a nearly quiet fluid's influence on g in the coupled problem is often lower-order in 
comparison to the other influences (including the primary influences driving the 
expansion). In other words, as long as the fluid behavior is tame, we expect that 
the coupled problem effectively behaves as a partially decoupled system in that g 
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can highly influence the fluid, but not the other way arouncQ. Thus, we have the 
following moral consequence of our Main Results: we expect that the addition of 
a perfect fluid (verifying p = c^p, < < 1/3) would not destroy the stability of 
a stable solution to a coupled Einstein-matter system, as long as the stable metrics 
are expanding at rates in the spirit of lAll - IA3I As we will soon discuss, this 
expectation has previously been confirmed when the expansion is generated by the 
inclusion of a positive cosmological constant in the Einstein equations; in this case 
the expansion is very rapid and verifies e^^*-* ^ e^*, where > is a constant. 

The key question from the previous paragraph (which is addressed by our main 
results) is whether or not the fluid behavior remains tame if it is initially tame. 
This is a highly non-trivial question whose answer depends on the expansion prop- 
erties of g. In the case of Minkowski spacetime (where e^*^*-* = 1), Christodoulou's 
monograph |Chr07| demonstrated that under a general physical equation of state, 
the constant solutions (i.e., p = p > 0, = Sq) to the relativistic Euler equations 
are unstabl(|l. More precisely, there exists an open family of initial data whose 
corresponding solutions develop shock singularities in finite time. Furthermore, the 
family contains arbitrarily small (non-zero) perturbations of the constant state, in- 
cluding irrotational data. We state some of Christodoulou's main results in more 
detail in Theorem 19.11 below. One very important aspect of his monograph is that 
it provides a complete description of the nature of the shock (we do not recall the 
fully detailed picture in Theorem 19. 1|) . His work can be therefore be viewed as 
a major extension of the well-known article }Sid85j . in which Sideris showed that 
singularities can form in the non-relativistic Euler equations under the polytropic 
equations of state for arbitrarily small perturbations of constant background solu- 
tions; in [Sid85] . a detailed picture of the singularity formation was not obtained 
because the argument involved the analysis of averaged quantities. 

In contrast to Christodoulou's shock formation result, the works |RS09j . |Spell| 
of Rodnianski and the author showed the following: when the scale factor verifies 
e^(*) ^ e^* for some "Hubble constant" H > 0, the equation of state isp = c^p, and 
< Cg < 1/3, then {g,p,u) (defined in (jl.2|) and p. 81) ) is a future-stable solution 
to the coupled Euler-Einstcin equations. The exponential expansion e^^*-*-* ^ e^^ of 
g in these solutions is driven by the inclusion of a positive cosmological constant A 
in the Einstein equationf0. where H = ^/A/3. These special "Friedmann-Lemaitre- 
Robcrtson- Walker" solutions {g,p,u) lie at the heart of many cosmological predic- 
tions. Furthermore, the case = 1/3 was recently addressed in jLVll] using an 
extension of Friedrich's conformal method (sec e.g. jFri02| ). which is closely related 
to the methods we use in Section [9l Roughly speaking, the conformal method is 
a collection of techniques for transforming the question of global-in-time existence 
into a question of local-in-time existence for the conformal field equations through 
changes of variables. This method often works when the energy-momentum tensor 
of the matter model is trace- free (as is the case of the fluid energy-momentum tensor 



This expectation is false for the Friedmann-Lcmaitrc- Robertson- Walker family of solutions to 
the Euler-Einstein system with no cosmological constant or additional matter model. In this case, 
the fluid itself is the primary influence driving the expansion, which is not accelerated. 

■^There was one exceptional equation of state which features small-data global existence when 
the fluid is irrotational; this equation of state leads (in the irrotational case) to the quasilinear 
wave equation for minimal surface graphs. 

^In the coupled problem, the scale factor corresponding to the explicit solution metric g verifies 
Friedmann's ODE, which leads to the asymptotic behavior e^'*' ~ e^*. 
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(|1.4p when = 1/3.) In summary, in the mild mitial condition regime, previous 
results show that exponential expansion suppresses the formation of fluid shocks 
when < < 1/3. 

We now mention another well-studied mechanism for generating accelerated ex- 
pansion in solutions to the coupled problem: the inclusion of a scalar field matter 
model in the Einstein-matter equations. In the simplest cases, the scalar field $ is 
postulated to verify a wave equation of the form {g"^)"^ DaDp^ = V'{^), where 
is a nonlinearity. As discussed in e.g. |Ren04a| . |Ren05a| . by adjusting V, one 
is able to exert control over the scale factor e^, at least for special solutions. Let us 
briefly discuss two interesting examples that lead to expansion rates verifying our 
hypotheses I All - 1 A3I First, in |Rin08| . Ringstrom studied nonlinearitics that verify 
V{0) > 0, y (0) = 0, y"(0) > O. For smaU $, the conditions on V generate a pos- 
itive cosmological-constant-like effect, and the coupled Einstein-scalar field system 
effectively behaves (for $ ~ 0) as if one had introduced the cosmological constant 
A = V^(O) into the Einstein equations. As discussed in the previous paragraph, this 
generates exponential expansion. The main resuh of |Rin08| was a proof of the fu- 
ture stability of a large family of background solutions to the coupled Einstein-scalar 
field system (with no fluid present). As a second example, in |Rin09] . Ringstrom 
studied the Einstein-scalar field system under the choice V{^) = Vbe~'^*, where 
Vq and A are positive constants. For certain explicit solutions, this choice of V 
generates power-law expansion e^^^*) ^ t'^ , where Q = -v/X, and the main result of 
[Rin09j was a proof of the future stability of these explicit solutions when Q > 1 
(i.e., when the expansion is accelerated). 

1.2. Outline of the analysis. 

1.2.1. The continuation principle. Our proof of the future stability of the explicit 
solutions ()1.8|) is based on a continuation principle (Proposition [33]) together with 
a standard bootstrap argument. With the aid of Sobolev embedding, the content 
of the continuation principle can be roughly summarized as follows: if the solution 
forms a singularity at time Tmax, then a certain order- iV (throughout the article, 

> 3 is a fixed integer) Sobolev norm Sn (see Definition 14. ip of the solution 
necessarily blows- up at time Tmax- Therefore, if one can derive a priori estimates 
guaranteeing that SN(t) remains finite for all t > 1, then one has demonstrated 
future-global existence. To show that SN(t) remains finite, we will assume that 

(1.9) SNit) < e 

holds on a time interval [l,r) of existence, where e is a sufficiently small positive 
number; by standard local existence theory, this assumption holds for short times 
when the data arc sufficiently small. Based on this assumption, we will then derive 
inequalities for S^it) leading to an improvement on the assumption (|1.9p . as long 
as the data and e are sufficiently small. Consequently, since S]\r{t) is continuous, it 
can never exceed a certain size, and future stability follows from the continuation 
principle. The complete argument is provided in the proof of Theorem l8.ll 

1.2.2. Energy currents. The only general method we know of for estimating so- 
lutions to quasilincar hyperbolic PDEs (such as the relativistic Eulcr equations) 
involves the derivation of coercive integral identities. This framework is often re- 
ferred to as the energy method, and its raison d'etre is connected to the following 
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basic difficulty: the time derivative of Sn (t) cannot be controlled in terms of the so- 
lution itself, for it is a crudely defined Sobolev norm which has no special structure 
tying it to the evolution of the fluid solution. To estimate S^it), we will construct 
coercive quantities SNit) called energies. SNit) will be used to control the norms 
of the up-to-top-order (i.e., up-to-A^*''— order) spatial derivatives of the solution via 
the divergence theorem (see (jl.lip ). We remark that the coercive nature of £'jv(i) 
is established in Proposition 16. II We also remark that in the case = 0, we will 
define separate energies EN-.veiodty and SN-i-.density for the four-velocity and the 
density. The reason is that in this case, the evolution of the four velocity decouples 
from that of the density, and furthermore, the density is necessarily one degree less 
differentiable than the four-velocity. 

To construct the energies, we will use a version of the method of vectorfield mul- 
tipliers. Over the last few decades, this fundamental method has been applied in 
many different contexts to a large variety of hyperbolic PDEs. All applications of 
this method are essentially extensions of Noether's theorem |Noe71| . In particular, 
wc make use of the multiplier method framework developed in |ChrOO| , a work that 
can be viewed as a geometric extension of Noether's theorem to handle regularly 
hyperbolic (see |ChrOO| for Christodoulou's definition of regular hyperbolicity) sys- 
tems/solutions that are derivable from a Lagrangian but perhaps lack symmetry 
(and hence the energy estimates contain error terms). The multiplier method allows 
us to construct compatible "energy" currents with the help of an auxiliary multi- 
plier vectorfield; we will further discuss the specific multiplier used in this article 
later in this section. The energy currents are special solution-dependent vectorfields 
J'^[ds'W,ds'W], and the energies are 



(1.10) £N{t)={ y / j°[dsW,dsW]d'^s 

Here, ds is an order < N spatial differential operator, W = (ln(e'^'-^+'^=^^p/ p), u^,u'^, u^) 
is an array of fluid variables, and the precise definition of J'^ is given in Definition 
14.21 As is suggested by the notation, J^[95W, cJ^W] depends quadratically on 
SgW. Furthermore, the coefficients of the quadratic terms depend on the solution 
W itself. 

In this article, we will not fully explain our construction of the currents J*^. 
Roughly speaking, these currents exist because the relativistic Euler equations are 
the Euler-Lagrange equations corresponding to a Lagrangian. For such PDEs, 
Christodoulou's aforementioned framework |ChrOO] explains how one can construct 
suitable currents J'^ using a version of the multiplier method. In fact, energy 
currents for the relativistic Euler equations were first derived by Christodoulou 
in the Lagrangian-coordinate framewor ifl in |ChrOO| . and later in the Eulerian- 
coordinate framework in |Chr07j . The Eulerian-coordinate fluid energy current 
framework has since been applied by the author and his collaborator in various 



contexts; see |Spe09a| , |Spe09b| , |Spell| , and [SSTl] 



^Interestingly, the relativistic Euler equations fall just outside of the scope of the methods of 
IChrOTI ; in the Lagrangian-coordinate framework, the methods of IChrOTI lead to semi-coercive 
estimates rather than fully coercive estimates. Nevertheless, we are able to derive fully coercive 
estimates in the Eulerian-coordinate framework. 
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Assuming that < < 1/3 and that (|1.9p holds for sufficiently small e, we will 
exploit the following two fundamental properties of : 

(1) j^[dsW,dsW]^ (9aln(e3(i+<=')^V//o) ) Ve^^^EjLil^aM^T (sec the proof 
of Proposition 16. ip . 

(2) For a solution W, it can be shown that 

df.iJ" [dsW, dsW]) = Fs{n{t),Lj{t);W, dW, d^^^W, 5^^^ W, • • • , a^'"'"^) W, ds 
The important point is that the smooth function does not depend on 
^(l"l+i)W^ where d denotes the spatial derivative gradient; see (|4.18p . 

These properties allow us to exploit the following basic integral identity, which 
is a version of the divergence theorem (see also ()4.22p and (|4.24ap - (|4.24bp ): 

(1.11) ^£Ut)= V / d^{j>^[dsW,d3W])d^x 

Above, Fn is a smooth function of its arguments. By (1), the right-hand side 
of (|l.lip is controllable in terms of £N{t)- Thus, £N{t) is a norm-like quantity 
whose time evolution can be estimated. Roughly speaking, in order to apply the 
continuation principle, we have to show that the identity (|l.lip prevents £N(t) 
and S^it) from blowing up in finite time. In particular, we have to carefully 
investigate the structure of Fn{- ■ ■); see the extended discussion in Section [1.2.41 
We also remark that in the case = 0, we will utilize two distinct currents Jyeiodty 
and Jl^^nsity ™ order to handle the partially decoupled nature of the relativistic 
Euler equations under the dust equation of state; this is of course connected to 
our previously mentioned use of separate energies £N-veiocity and £N-i-density when 

It is interesting to note that Christodoulou showed |ChrOO| that the Lagrangian 
for relativistic fluid mechanics is generic, which means that up to null currentf0 
(which by themselves lead to non-coercive integral identities) the method of vector- 
field multipliers generates all possible integral identities that a solution can verify. 
Furthermore, he showed that a certain widely used subclass of compatible cur- 
rents, namely those generated by domain multiplier vectorfields, have the following 
property: the only useful members of this class (i.e., generating an L^— coercive 
compatible "energy" current) are generated by domain multiplier vectorficlds that 
belong to the inner core of the characteristic subset of the tangent space. Now for 
the relativistic Euler equations, the inner core is the degenerate one-dimensional 
set span{u} (see |Spe09a| ); this is connected to the fact that the fluid vorticity 
two-form verifies a transport equation in the direction of u. Thus, it can be shown 
that up to multiplication by a scalar function, the only coercive compatible energy 
current generated by the domain multiplier vectorficlds is the one defined below in 
Definition\4-.2\ 



null current is a solution-dependent vectorfiold whose divergence can be shown to be in- 
dependent of the solution's derivatives without having to use the equations of motion (which are 
the relativistic Euler equations in the present article) . Such null currents have a special algebraic 
structure and have been completely classified |ChrOO| . 
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1.2.3. The analytic effect of expansion and the restriction < < 1/3. From the 
analytic standpoint, the stabihzing effect of spacetime expansion can be summarized 
as follows: the scale factor e^^*^ generates the dissipative w(i)(3Cj — 2)u-' term on the 
right-hand side of p.7cp . The sources of these dissipative terms are the Christoffel 
symbols of g relative to the coordinate system {t ), which are provided in 
Lemma 15751 In the cases < < 1/3, equation p.7cp can be rewritten as 

(1.12) dt[e^.'^{n)u^] = u{?,cl - 1 + :JM|[e^^^(r!)uJ] +e^^^(17)A'^ 

where ^ is the function from the hypotheses IA2I - IA3I and u;(i) ■^il{t). If we 
assume that the right-hand side of (|1.12p is negligible (we will in fact establish 
this assumption when < < 1/3 and (|1.9p holds), then by integrating the 
"ODE" (fLT2l) . we deduce that \u^\ < ee~^^^*^ {n{t)) . Therefore, the size of the 
spatial part of u as measured by g is \gabu"'u^\ = e^^5abU°'u^ ^ {^{t)) . Using 

hvpotheses lAll - [A3[ it follows that when < < 1/3, the quantity gabu'^u^, which 
appears e.g. as an inhomogeneous term in the Euler equation p.2a|) . is expected 
to decay towards its background solution value of 0. Furthermore, (see (jl.6p ) 
is expected to decays towards its background value of 1. In the case = 0, our 
assumption (|1.9p will directly imply that \gabu'^u^\ ^ e^e^^^(*\ and we again expect 
decay towards 0. In the case = 1/3 (in which ^ = 1), \gabu°'u''\ is expected to 
remain at size e^. In contrast, if > 1/3, this heuristic analysis suggests that 
there is nothing preventing the quantity \gabu'^u^\ from growing unabatcdly; we 
thus anticipate that there may be instability in these cases. We remark that in 
|Ren04bj , Rendall also detected evidence of instability for > 1/3 through the use 
of formal power series expansions. 

1.2.4. Comments on the analysis of (|l.lip . Much of our work goes into estimating 
the J^,FN{n{t),uj{t);W,dW,d'-^^W,d!'^^W,--- ,d!'^^W)d^x term on the right- 
hand side of (|l.lip in terms of the Sobolev norm Sn under the bootstrap assump- 
tion (|1.9p . The end result of these estimates in the cases < < 1/3 is inequality 
(|5.10p . and in the case = is the inequalities (|5.8p - (|5.9p . Our analysis is 
based on the identity (|4.18p for 9^ J'' (in the cases < < 1/3), the identities 
(|4.19ap - (|4.19bp for <9^i j^eZocity ^.nd d^jj^^^^^^^ (in the case =0), on standard 
Sobolev-Moser type estimates (see the Appendix), and on Lemma |4?2] Lemma l472l 
is a particularly important ingredient in this analysis, for this lemma dissects the 
algebraic structure of the inhomogeneous terms in the equations verified by the 
quantities d^W. Understanding this structure is essential since these inhomoge- 
neous terms appear in the expression (|4.18p for 9p( J^[95W, SsW]) (and hence in 
F/v(- • • )), and also in the analogous expressions (|4.19ap - (|4.19bp in the case = 0. 
We remark that the quantities dsVi for |a| < verify linear (in dsW) PDEs 
whose principal coefficients depend on the solution W and whose inhomogeneous 
terms depend on the derivatives for \(3\ < \d\. These PDEs are known as the 
equations of variation; we describe them in detail in Section [4.21 

We now discuss some of the subtleties in our derivation of the estimates (|5.8p 
- (|5.10p which are the main estimates used in our derivation of the key energy 
inequalities (|7.1ap - (|7.3p of ProDOsition l7.1l For simplicity, we will focus the present 
discussion on the cases < < 1/3 and the estimate (|5.10p . which provides an 
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upper bound for J2\a\<N /k^ '^tJ-i^'^i^^s^ , ds^]) <Px in terms of the norm Sn- A 
very important aspect of our argument in these cases is the following: in order 
to derive (|5.10p and to close our global existence argument, we need to prove that 
the lower-order spatial derivatives of decay faster than its top-order derivatives. 
To see why this is the case, we will now show that without this improved decay, 
there would be an obstruction to global existence. We illustrate this obstruction by 
closely investigating a specific term from the right-hand side of (|l.lip . In particular, 

the -jrj-j^lyL term from the identity (|4.18|) below, where ^ = —uj ffab""" is the 
term on the right-hand side of p.2ap . leads to the identity 

d r 1 

(1.13) -El{t) = -2cl UJ— ^ u''u''Ld^x + --- , 

at Jg_3 u 

where L = In (e3(i+-?)f^p/p) , and S^^it) « + e^^^ ^^^i Using 

(|L9l) . (fTTSl) . the Sobolev embedding resuh \\u^\\l'=- < \\u^\\h^ < e^^^fw, and the 
estimate u° > 1, we deduce that 

(1-14) ^^£Ut)<u^£Ut) + -- - ■ 

Inequality (|1.14p allows for the dangerous possibility that Spfit) grows unabatedly, 
and is therefore not sufficient to prove our main stability theorem. 

dcf 

To remedy this difficulty, we introduce the lower-order Sobolev norm Un-i = 

/ 3 \l/2 

e^(f2)( J2j=i ) (see Definition l4.1| ). Here, ^ is the function from the 

hvpotheses I A2I - 1 A3I The main point is the following: if we knew that Un-i < e, 
then the Sobolev embedding estimate ||u-'||loc < Hu-'Hh^ ^ e~^I^~^{^)UN-i ^ 
ee~^^^~^(r2) would allow us to upgrade (|1.14p to 

(1.15) j^EUt) < euj{t).^-\m)£Ut) + ■■■■ 

Because of our assumptions ! All - IA3l on the scale factor, we have that ijj{t)sS^^{Vl{t)) 
= ^1 dfij e L\{[l,oo)). Therefore, the first term on the right-hand 

side of (|1.15l) is amenable to Gronwall's inequality. Assuming that the remaining 
terms • • • could be treated similarly, we would therefore be able to derive an a priori 
estimate for iS^(t) guaranteeing that it remains uniformly small for all time; thanks 
to the continuation principle, this is the main step in proving global existence. At 
its core, the proof of our main stability theorem is essentially a more elaborate 
version of the estimate (jl.lSp . 

It remains to discuss how we obtain the desired control over the norm Un~i- 
Roughly speaking, we will use the improved decay for the suggested by equation 
(|1.12p . It turns out that the error terms corresponding to this equation are small 
enough that we can effectively control the time derivative of iYjv-i(i) even if we 
only have control over the weaker norm '^\s\=n 6^II^5^"' IIl^ of the top-order 

derivatives (which will follow from the bootstrap assumption (|1.9p ). The precise 
estimates corresponding to this fact are provided in inequalities (|4.23p and (j5.4dp 
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below. We emphasize that in order to derive the improved decay (improved rela- 
tive to the rate predicted by the Sobolev norm bootstrap assumption (|1.9[) ). what 
matters is the interaction between the hypotheses lAll - IA3I on the scale factor 
and the structure of the error terms. Note that (|1.12p is not literally an ODE in 
since spatial derivatives of the solution are present on the right-hand side as 
"error terms." Therefore, the approach we have described only allows us to derive 
improved decay estimates for the below-top-order derivatives of the solution. We 
remark that wc do not use the norm Um-i in our analysis of the cases ~ 0, 1/3. 

1.2.5. Sharp results in case ~ 1/3. Our sharp analysis in the case — 1/3 is 
based on the well-known conformal invariance of the relativistic Euler equations 
in this case; see Section [9] for more details. Roughly speaking, when = 1/3, one 
can perform a changc-of-time- variable ^ = e~^*^*\ t(1) = 1, and also a rescaling 
of the fluid variables in order to translate the problem of interest into an equivalent 
problem in Minkowski spacetime. Now if assumption I A 1 1 holds but e"^^*-* ds = 
oo, then the analysis of Section |9] shows that the instability of the background 
fluid solutions follows easily from the instability results derived in |Chr07| : see 
Corollary 19.21 Furthermore, the arguments given in Section [9] could easily have 
been extended to show the future stability of the background fluid solutions when 
Cj = 1/3, hypothesis |Al] holds, and /J^^ e"^^^-* rfs < oo. However, wc instead 
chose to use the energy-method/continuation principle framework to prove future 
stability (see Theorem 18. ip since this framework is stable and since there is no 
known alternative to this framework in the cases cj. G [0, 1/3). 

1.3. Outline of the article. The remainder of the article can be summarized as 
follows. 

• In Section [2] we introduce some standard notation and conventions. 

• In Section [3l we derive an equivalent version of the relativistic Euler equa- 
tions that is useful for our ensuing analysis. We also introduce some stan- 
dard PDE matrix-vector notation. 

• In Section m we introduce the Sobolev norms and the related energies that 
we use to analyze solutions. 

• In Section [SJ we derive the Sobolev estimates that play a major role in our 
derivation of differential inequalities for the fluid norms and energies. 

• In Section [6] we prove a simple comparison proposition, which shows that 
the energies we have defined can be used to control the norms. 

• In Section [3 wc use the estimates of Section [5] to derive differential inequal- 
ities for the norms and energies. 

• In Section [8l we use the differential inequalities to prove our main future 
stability theorem. 

• In Section ^ we use a well-known result of Christodoulou to show that in 
the case = 1/3, the non-integrability condition e~^^^^ ds = oo leads 
to the nonlinear instability of the explicit fluid solutions. 

2. Notation 

2.1. Index conventions. Greek "spacetime" indices a,/3, ••• take on the values 
0,1,2,3, while Latin "spatial" indices a, &, • • • take on the values 1,2,3. Pairs of 
repeated indices arc summed over their respective ranges. We lower spacetime 
indices with g^i, and raise them with {g~^)^'^ ■ 
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2.2. Coordinate systems and diflferential operators. We perform most of our 
computations relative to the standard rectangular coordinate system (x", ,x^,x^) 
on M.^. We often write t instead of x^. In this coordinate system, the metric g is of 
the form (|1.2p . We often use the symbol 9^ to denote the coordinate derivative ^|jr , 
and we often write dt instead of Oq. Occasionally, (sec e.g. the derivation of (j5.27p ) 
we also use the rcscalcd time variable t introduced in Section [HI The corresponding 
rescaled partial time derivative is denoted by dr- All spatial derivatives are always 
taken with respect to the coordinates {x^,x'^,x^). 

If a = (ni, 712, na) is a triplet of non-negative integers, then we define the spatial 

multi-index coordinate differential operator ds by ds = ^"^Sj^Sg^. We use the 
notation |a| ni + 112 + to denote the order of a. 



(2.1) 

r s 

a— 1 a— 1 

denotes the components of the covariant derivative of a spacetime tensorfield T. 
The Christoffel symbols T"^^ are defined in p. lip . 

^(Af) jwi - iy^ denotes the array containing of all of the iV*'' order spacetime coordi- 
nate derivatives (including time derivatives) of the component T"^;^.'.'.'^^ ■ ^'■^■'T^^.'.'.'^^ 
denotes the array of all TV*'* order spatial coordinate derivatives of the component 
r^^.;;^^. We omit the superscript when A'' = 1. 

2.3. Norms. We define the standard Sobolcv norm of a function ||/|| „ as follows: 



(2.2) ||/L«=YE / \ds.fit,x\x^x')\'d'x 



S\<N' 



1/2 



The above volume form d'^x corresponds to the standard flat metric on MP. 
Wc denote the N*^ order homogeneous Sobolev norm of / by 



(2.3) lls^'^VlL. = E ll^-s/IL- 

\S\ = N 

If C M" or C T", then (J^) denotes the set of iV-times continuously 
differentiable functions (either scalar or array- valued, depending on context) on 
the interior of M. with bounded derivatives up to order N that extend continuously 
to the closure of We define the norm corresponding to this function space by 



(2.4) |FU,^'=' 5^ esssup|a,^F(.)|, 

/|<Af 

where df is a multi-indexed differential operator representing repeated partial dif- 
ferentiation with respect to the arguments . of F, which may be either spacetime 
coordinates or metric/fluid components depending on context. When A'^ = 0, we 
use the slightly more streamlined notation 
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(2.5) |FU'='esssup|F(.)|. 
Furthermore, we define 

(2.6) |F(^)L,^=^^ 

\I\=N 

When .8 = R'^, we sometimes use the more familiar notation 

(2.7) ||F|U^'='esssup |F(x)|, 

(2.8) ll^llcf = E \\9»F\\l^- 

\a\<N 

If A is an m X n (often 4 x 4 in this article) array- valued function with entries 
Ajk, (1 < J < 1 < fc < n), then in Section[Sl we write e.g. |jy4.|jjLfN to denote the 
mx n array whose entries are ||Ajfe||/jjv. We use similar notation for other norms 
of A. 

If / C M is an interval and X is a normcd function space, then (I, X) denotes 
the set of iV-timcs continuously diffcrentiablc maps from / into X. 

2.4. Numerical constants. C denotes a numerical constant that is free to vary 
from line to line. Wc sometimes write e.g. C{N) when we want to explicitly indicate 
the dependence of C on quantities. We use symbols such as c, C,, etc., to denote 
constants that play a distinguished role in the discussion. Wc write X <Y when 
there exists a constant C > such that X < CY. We write X Ri Y when X < Y 
and Y <X. 

3. Alternative Formulation of the Relativistic Euler Equations 

In this section, we derive an equivalent version the relativistic Euler equations 
(jl.7ap - (|1.7b[) under the equation of state p = c^p; wc will work with this version 
for most of the remainder of the article. We then briefly discuss a classical local 
existence result and a continuation principle. The latter provides sufficient criteria 
for the solution to avoid forming a singularity. 

3.1. Alternative formulation of the relativistic Euler equations. Our al- 
ternative formulation of the relativistic Euler equations is captured it the next 
proposition. 

Proposition 3.1 (Alternative formulation of the Euler equations). Assume 
that the fluid verifies the equation of state p ~ c^p. Let 

(3.1) L*^ln(e3(i+^')^V/p) 

be a normalized energy density variable, where p > is the constant corresponding 
to the background density variable p~ pe"'^'^^'^^'^^ Then in the coordinate system 
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(<, x^, cc^, x"^), the relativistic Euler equations (jl.7ap - (jl.Tbp are equivalent to the 
following system of equations in the unknowns (i,u^,u^,u^) (for j = l,2,3j; 



(3.2a) + (1 + cl) (l-)uadtu'' + (1 + c^dau'^ = -u^^^^^gabu'^u' 



(3.2b) + '—W^daL = w(3c2 - 2)u'^u^ . 

(1 + ci) 

Above, 



u° = (1 + 5afc"°"')'/', = ^^^^^^ + (.9"')^^ ^(t) = -f^(i). 

Furthermore, u'^ is a solution to the following equation: 



(3.3) + -— ^nO°9„L = w(3c2 ~ l)gabu''M^ 

Proof. To obtain p.2ap . we first expand the covariant differentiation in (|1.7a[) to 
deduce the following equation: 

(3.4) u"9„lnp+(l + c2)a„ii" = -il + cl)T^"^u^. 

Lemma 13.31 implies that ^^pU^ = 3uju^, while equation (jl.6p implies that dtu'^ ~ 
:^{uadtu°' + wgafcu'^u''} . Equation p.2ap now follows from these identities and the 
identity 9t In p = 9t In L — 3(1 + c'^)uj. 

Similarly, to obtain ()3.2bp , we first expand the covariant differentiation in (jl.7bp 
to deduce the following equation: 



(3.5) u"d^u^ + —^W-d^ \np= -rj' 

(1 + ci) « P 

2 

Equation (|3.2b|) now follows from p.5p . Lemma [?751 and the identity jj^f^Tl-'°'da Inp 

2 

■jj^^W^daL — Sc^wu^u^ . The proof of (|3.3p is similar, and we omit the details. 

□ 

Remark 3.1. Note that in terms of the variables {L,u^ ,u'^ ,u^), the background 
solution (11.81) takes the form 



(3.6) (L,u^u^u^) = (0,0,0,0). 

Part of our analysis involves solving for dtL and dtu^^ and treating spatial deriva- 
tives as inhomogeneous error terms. This approach can also be used to analyze dsL 
and daU^ 1 for \ci\ < N — 1. As a preliminary step in this analysis, we solve for dtL 
and dtu^ ■ 



The Stabilizing Effect of Spacetime Expansion 



15 



Corollary 3.2 (Isolated Time Derivatives). Let {L,u^ ,u'^ ,u^) be a solution to 
the relativistic Euler equations (|3.2ap - (j3.2b[) , where is defined by (jl.6p . Then 
the time derivatives of the fluid variables can be expressed as follows: 

(3.7a) dtL = A', 

(3.7b) dtu" = A'", 

(3.7c) dtu^ = uj{3cl - 2)u^ + A'^ 

where the error terms A', A'", and A'J are defined by 



(3.8a) 



2\n q„2nJ 



1 - c 



s (tiO)2 

2 9abU''u''y\ 2 2N^aU" , (1 + C?) 



(3.8b) 



A'O = uj{3cl - l)\ , \ - ^^^\ Z^}^ ^u-daL 



„2 gab'u'^u'' ^ gabu'^u''dku'' 



2 _ ^2 gabW^U" I -| _ „2 gabW^u" 

(3.8c) 



gabu°'u . ^ dgu'^ gabU^u dh'U 



A'^- = a.c^(3c^ - l)un ^^^^ + c^-^- ^^^^ ^ 



r B „J 

[u ) 

Proof. The proof consists of tedious but simple calculations; we omit the details. □ 

Remark 3.2. Note that certain terms on the right-hand side of p.8ap - (|3.8cP 
vanish when = or = 1/3. This vanishing is important in our analysis of 
these cases. 

In the next lemma, we provide the ChristofFel symbols of the metric g relative 
to our coordinate system. The lemma was used in the proof of Proposition 13. II 

Lemma 3.3. The non-zero Christoffel symbols of the metric 
g = -dt^ + e2f^(*) Y.]=i{dx'f are 

(3.9) r/, = r,°^. T^o = T^\^uj51 (j-fc = l,2,3), 

where 



(3.10) 
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Proof. The lemma follows via simple computations from the definition 



(3.11) -(^g-^)^\d^gxu + d^g^x - d>.g^,). 

□ 

3.2. Local existence and the continuation principle. In this section, we state 
a standard local existence result for the system p.2ap - (j3.2bp . 

Theorem 3.1 (Local Existence). Let N > 3 be an integer. Let L = i|t=i = 
ln(p/p)|j^i, — u^\t=i, {j = 1,2,3), be initial data for the relativistic Euler 
equations (j3.2ap - (j3.2bp satisfying 



(3.12a) L e 

(3.12b) e H 



H^-\ c2 = 0, 
H^, < c2 < 1/3, 

-N 



where p > is a constant. Assume that sup^ggs \L\ < oo. Then there exists a real 
numberT^ > 1 such that these data launch a unique classical solution (L, , , u^) 
to the relativistic Euler equations (|3.2ap - (j3.2bp existing on the spacetime slab 
[1,T~^) X M^. Relative to the coordinate system {t ), the solution has the 

following regularity properties: 



(3.13a) 
(3.13b) 



= r c"([i,r+)xR3), 
\ ci([i,r+)xR3), 

e C^[l,T+) X M3), 



c2 = 0, 

< c2 < 1/3, 



(3.14a) L e 



^([1,^+),^^-!), cl^O, 
C°{[1,T+),H^), 0<c2<l/3, 



(3.14b) w°-l, w^' e C"([l,T+),iJ^), 

In addition, there exists an open neighborhood O of {L,u^) such that all data 
belonging to O launch solutions that also exist on the slab [1,T_|-) X R3 and that 
have the same regularity properties as (L,u^). Furthermore, on O, the map data ^ 
solution is continuous. 



Proof. Theorem 13.11 can be proved using a standard iteration or contraction map- 
ping argument based on energy estimates for linearized equations that are in the 
spirit of the estimates derived in Scctions|6]and[7l See e.g. jHor97| Ch. VI], |Spe09b| 
for details on how to use such energy estimates to deduce local existence. □ 



In our proof of Thcorcm lS.ll wc invoke the following continuation principle, which 
provides standard criteria that are sufficient to ensure that a solution to the rela- 
tivistic Euler equations exists globally in time. 

Proposition 3.4 (Continuation Principle). Assume the hypotheses of Theorem 
\3.1l Let Tmax be the supremum over all times such that the solution (L, u^) 
exists on the interval [1,Tj^) and has the properties stated in the conclusions of 
Theorem ] 3. 11 Then in the case = 0, if Tmax < oo, we have 
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(3.15) 



lim sup |||L(s,.)|lL~+^(|k^"(s,-)| 



1 > = oo. 



Furthermore, in the cases < Cg < 1/3, if T^ax < oo, we have 



(3.16) 



lim sup ■ 

t^T;;,^^ o<s<t 



\L{s, 



3 

E 



(ii«^'(5,-)ii( 



oo. 



Proof. See e.g. [Hor97j . |Spe09a| for the ideas behind a proof. The case = is 
special because in this case, the evolution of the decouples from that of L, and 
furthermore, the evolution equation p.2ap is linear in L. □ 



3.3. Matrix-vector notation. It is convenient to write the system p.2ap - ()3.2bp 

in abbreviated form using standard PDE matrix- vector notation: 



(3.17) 



Here, W is a fluid variable column array, and b is a column array of inhomogeneous 
terms: 



(3.18) 



W 



u 


( 


, dcf 

b = 




V 



a;(3c2 - 2)u°u^ 
w(3c2 - 2)u°u'^ 



The are 4x4 matrices defined by 



(3.19a) 



(3.19b) 



( u° 




(n 




■\ til 


(1 




+ 






nio 


















1(1+^^); 


n2o 


















/ <?. ^ 


n3o 




















(1 


+ C' 















TTll 

(1+c^)^^ 


















c„ tt21 
(1+c^) 







ui 











V(l+c=)n 



















and analogously for A^, A^. Furthermore, for later use, we calculate that dct(^°) 



(.")2{(uY(l-cD + c^}, 
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(3.20) 



oo\-i 



n 



20 



30 



1*0 -di 



and di 



dcf 



(1 + 


„2\ M2 


(1 + 


„2 \ M3 \ 




ni0u2 












-d2 




n20u3 


(tlU)2 


cl 


n30u2 




-^3 / 



4. Norms, Energies, and the Equations of Variation 

In this section, wc define some Sobolev norms and related energies, all of which 
arc used in the proof of our main stability theorem. Our choice of norms is mo- 
tivated in part by the continuation principle (Proposition 13. 4| : with the help of 
Sobolev embedding, our norms will control the quantities appearing in the state- 
ment of the continuation principle. The energies are introduced in order to control 
the up-to-top-order spatial derivatives of the solution. Unlike the Sobolev norm Smi 
the energy's time derivative can be controlled in terms of the energy itself with the 
help of the divergence theorem. This is because our energy is defined with the help 
of energy currents J'^, which are solution-dependent, coercive vectorfields whose 
divergence can be controlled; see the discussion in Section ri.2.21 We also introduce 
the equations of variation, which are the PDEs verified by the spatial derivatives 
of the solution. The structure of the equations of variation plays an important role 
in our derivation of estimates for the spatial derivatives. 



4.1. Norms. 



3\T 



Definition 4.1 (Norms). Let iV be a positive integer, let W = {L , , v? , v!^) 
be the array of fluid variables. In the cases < < 1/3, we define the lower-order 
fluid velocity norm lAi^^i{t) > by 



(4.1) 

3 3 

Z^Ar_/='e^^^(17)(^||u^"||^„_,)'^'-e"^(17)(^ ^ / {d^u^fd'xf'^ 

J = l 3 = l\S\<N-l''^^ 

In the cases < < 1/3, we define the fluid norm SN{t) > by 



\L\\H--.+e^'^j:' \\u^\h-, 



(4.2) Sn ={ \\L\\hn +e''j:'\\u^\HN +Un-i, < < 1/3, 

\L\\H-+e''Ej=ih'\\H-, c2 = l/3. 

Above, ^ is the function from the hvpotheses I A2I - 1 A3I of Section [TJ 

In the case = 0, we also define the fluid velocity norm SN;veiocity{i) > by 



(4.3) 



*-^N;velocity — ^ 
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Remark 4.1. The proof of our main stability theorem will show that all of the 
norms remain uniformly small for all future times if they are initially small. Fur- 
thermore, in the cases < < 1/3, the boundedness of Um-i shows that the 
lower-order derivatives of have an norm that decays at least as fast as 
e~^,^~^(ri), while the boundedness of Sm shows only that the top-order deriva- 
tives (i.e. \a\ = N) of decay at least as fast as e~^^; sec Section [1.2.41 for a 
discussion of the importance of the improved decay for the lower-order derivatives 
of the . We remark that in the special case = 0, the boundedness of Sn shows 
that all spatial derivatives of u> decay in at least as fast as e~^^. Note also that 
in the case = 1/3, we do not prove an improved decay rate for the lower-order 
derivatives of (and hence the IAm-i norm is not needed for the analysis in this 
case). Our inability to show improved decay in this case is intimately connected 
to the fact that when = 1/3, the — (Sc^ — 2)uju^ term on the right-hand side of 
equation (j3.7cp suggests that we can only prove that the lower-order derivatives of 
decay like (i.e., the same decay as for the top-order derivatives). 

4.2. The equations of variation and energies for the fluid variables. In this 
section, we define the fluid energy. This energy, which depends on the up-to-top- 
order spatial derivatives {dsL, d^u^ ,ds_v?, dsu^) \d\ < N, plays a central role in the 
proof of our main stability theorem. We note the following previously mentioned 
exception to the preceding sentence: when = we will only be able to control 
dsL for |c5| < — 1. In order to study the evolution of the solution's derivatives, 
we will have to commute the equations (|3.2ap - (j3.2b[) (or equivalently equation 
(jS.lTp ) with the operator ds- The quantities {dsL.dsu^ ,dsv? ,dav?) verify hnear 
(in (9(5£, dsu^ jdsu^ ,dsu^)) PDEs with principal coefficients that depend on the so- 
lution and inhomogeneous terms that depend on {d^L, d^v} ,d^v? , d^u^), < 
see Lemma [321 for the details. We refer to this system of PDEs as the equations of 
variation, while the unknowns {L^v} ,u^) '= {dsL,dsu^,dsu^,dsu^) are called 
the variations. More specifically, we define the equations of variation in the un- 
knowns (L, u^, u^, ?i'^) corresponding to (L, u^, w^, u^) as follows: 

Equations of Variation 



(4.4a) u''dc.L + {l + c^^)(^^yadtu'' + {l + cl)dau'' 

(4.4b) + '—WdaL = ujiicl - 2)w"w^ + & . 

(1 + cf) 

The terms ^, a;(3Cs — 2)u^u^ , and 0^ denote the inhomogeneous terms that arise 
from commuting p.2ap - p.2bp with ds- Note that we have split the inhomogeneous 
term in (j4.4bp into two parts; the a;(3c^ — 2)u^u^ term is primarily responsible for 
creating decay in ii^ , while & will be shown to be an error term. 
Using matrix-vector notation, we can rewrite (|4.4ap - (j4.4b[) as 



(4.5) A^^dpW = I, 

where them matrices arc defined in (|3.19ap - (j3.19bp . and 
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(4.6) 

(4.7) I Lo{3cl - 2)^0(0, w\ u^y^ ^ gi^ (g2^ 

To each variation (u^, u^, it'^), wc associate a quantity u*^ defined by 



(4.8) 



. n dcf i n ■ b 



This quantity appears below in the expression (|4.16p , which defines our fluid energy 
current. The importance of the definition (|4.8p is that it leads to the identity 
gapu°'u^ = 0, which is essential for the derivation of the divergence identity (|4.18p 
below. In our analysis, we will need the following lemma, which essentially states 
that V? is a solution to a linearization of p.3p around (L, u^, it^, u"^). 

Lemma 4.1. Assume that {L,u^ ,u^,u^) is a solution to the equations of variation 
(|4.4ap - (j4.4bp corresponding to [L,u^ ^u^ ^u^), and let iiP = -^Uau"" be as defined 
in (|4.8p . Then vP verifies the following equation: 



(4.9) 

where 

(4.10) ©° = gab 



(l + c2) 



Proof. By the definition of vP, the left-hand side of (|4.9p is equal to 



(4.11) + --^^n°"9„i + 



■if + liogabU^u 



On the other hand, contracting equation (j4.4bp against Uj and using the identity 
u„n"" = -iion°" = M°n°", we conclude that 



(4.12) UaU^'do^il'' + ^^n0"9,i = (3c2 - 2)c^5„fcU''u'' + ffafaU"©'. 

Multiplying (|4.12p by ^^rr and using (|4.1ip . we arrive at (g^. □ 

In the next lemma, we investigate the structure of the inhomogeneous terms in 
the equations of variation verified by a solution's derivatives {d^L, dsu^, dsu^, dgu^)^ ■ 
We again split the inhomogeneous terms into a decay-inducing piece and a small 
error term b^g. 

Lemma 4.2. LetW = {L tU^^vP^u^)'^ he a solution to the relativistic Euler equa- 
tions ([XTT]) . I.e., Al^dpW = b, 

b "Jcoi^-^^gabU-u', (3c2-2)«%i, {3c^^-2)u"u^ {3cy2)u°u'Y . Then (L, u\u^,u^f 

= {dsL,dsu^ ,dsu^,dsu^)^ is a solution to the equations of variation (|4.5p with 
an inhomogeneous term I that can be expressed as follows: 
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(4.13) I-bs + bAs, 

where 



(4.14a) hs = w(3c2-2)u°(o,95m\55u^95u^)^, 

(4.14b) bA5 {ds. ^h, ©I, ©1)"^ = {dsh - bs} + [A'^dsiiA^r'h] - ^^b} 

+ A"{(A°)-i^'^aaa5W - dsiiA'^r'A^'daW]}. 

Furthermore, iiP ■^^Ua'ii'^ is a solution to equation (j4.9p wit/i an inhomogeneous 
term 25ji defined by 

(4.15) 

froo/. Equations (|4?T3)) - (j4.14bp are a straightforward decomposition of the inho- 
mogeneous term I = A^da{{A^)-^Af'd^MV] + A%A^)-^ A^'^^,, da\W 
= A^ds{{,A^^y^h] + A°[(A")-iA°aa,5s]W, where [•, •] denotes the commutator. 
The relation (|4.15p follows directly from (|4.10p . 

□ 

4.2.1. The fluid energy currents. Our energy, which will control the up-to-top-order 
spatial derivatives of the solution, will be defined with the help of energy current 
vectorfields J^. 

Definition 4.2 (Currents). In the cases < < 1/3, to each variation W — 
(i, M^, u^, u"^)"^, we associate the following energy current, where ^"au" : 

(4.16) j" = j^^L^ + 2c2ij^i + (1 + ci)u^ga.pu"u^. 

(1 + cf) 

In the case = 0, we associate both a velocity energy current and a density 
energy current: 

(4.17b) j',eusUy\L,L]'^u^L\ 



Similar currents have been used in |Chr07j , |Spe09a| , Spe09b , and |Spell| . See 
the discussion in Section [I.2.2l 

Remark 4.2. We use the notation J'^[.,.] when we want to emphasize that 
depends quadratically on the variations. 
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4.2.2. The coordinate divergence of the currents. In this section, we provide expres- 
sions for the coordinate divergence of the currents. These expressions wih play a key 
role in our analysis of the evolution of the energies £n, £N;veiocity, and £N-i;density 

Lemma 4.3 (Divergence of the Currents). Let W = {L,u^ ,u^)^ be a 
solution to the equations of variation (|4.4aP - (j4.4b[) . Then in the cases < < 
1/3, the coordinate divergence of can be expressed as follows: 



(4.18) a^(>[W, W]) = j^^^L' + (1 + cl){d,u'^)[- {iiy + gam'^u" 



u L + Ac ui — L 

J ' 

+ 2(l + c2)(3c^-l)cjg,,wV 



u 



<Oifcl< 1/3 



(1 



L - 2(1 + ci)(S''u" + 2(1 + ci)gab'S''u'' 



In the case c^ = 0, the coordinate divergences of Jy^i^^iiy and J'Ji^nsUy '-^'^ 
respectively expressed as 

(4.19a) 

+ 2we2f^(w° - l)gabu''y' - 2ue^'\%uy 
- 2e2"©°w° + 2e'''gab^^u\ 
(4.19b) 54 j,^,„,,,^[i, L]) = {d,u^)L^ + ^^ug^bW^u^L 

^ gabu'^'S'L + -l-g,,u^{u''dku'')L 



(w0)2^"" (w0)2- 

- 2{daU°')L + 2^L. 

Remark 4.3. We stress that the right-hand sides of (|ITT5)) and (|4.19ap - (|4.19b|) 
do not depend on the derivatives of the variations. This property is essential for 
our analysis of the energies. We have organized the terms on the right-hand sides 
of and (|4.19ap - (|4.19bp in order to help us deduce the estimates ([O)) - ([CT)) 

and (jSTOl . 

Proof. To deduce (I4.18p . one can take the divergence of (|4.16p and then use the 
equations (|4.4ap - (j4.4bp and (|4.9p to replace the derivatives of W with inhomoge- 
neous terms. We leave the tedious details up to the reader. The proofs of (|4.19ap 
- (|4.19bp are similar but simpler. □ 

4.2.3. The definition of the fluid energies. We are now ready to define our fluid en- 
ergies, which will be used to control all spatial derivatives of the solution, including 
those of top order. 

Definition 4.3 (Energies). Let W (L, u^, u^, it'^)^ be the array of fluid vari- 
ables, and let >[•,•], j'^^i^^uyl- jdens^ty[■r] be the currents defined in gH]), 
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(|4.17ap - (|4.17bp . In the cases < < 1/3, we define the fluid energy £N{t) > 
(see inequahty (|6.4I) ) by 

(4.20) £N{t)= V / j'^[dsW,dsW]d'^x. 

In the case = 0, we define the energies £n -velocity (t) > and £N-i;density(t) > 
Oby 

(4.21a) £kueiocUy{t)= E / j°eiocUy[ds{^\ , ds{u\ u')] X , 

|a| < A' 

(4.21b) ^N^Udensityit) ~ E / "^densityi^sL , dgL] d^ X . 

|a|<A' — 1 

In the next corollary, we provide a preliminary estimate for the time derivatives 
of some of the fluid-controlling quantities. These estimates are the starting point 
for our proof of Proposition 17.11 which provides the integral inequalities that form 
the crux of our future stability theorem. 

Corollary 4.4 (Preliminary Norm and Energy Inequalities). Let W 

{L,u^,u'^,u^)'^ be a classical solution to the relativistic Euler equations p.l7p . Let 
U]\[-i(t)^ 5jv(i), £N:veiocity, o,nd £ N -1-density be the norm and energies defined in 
Definitions \4-l\ and \4-3\ Then in the cases < cj. < 1/3, the following identity 
holds: 



(4.22) tA^n)= E / d,{j'^[d3W,dsW])d'x. 

In the cases < cj. < 1/3, the following differential inequality holds: 
(4.23) 

j^iUl^,) < 2{3cl - 1 + ^}coUj,_, + E e"^(f^)|| A-||^»-.. 

In the case = 0, the following identities hold: 

(4.24a) f^{£l,,,eloc^ty) = E [ ^ d^{j'Vsi^\ u^) , dsiu\ IL^)]) d^ X , 
d f 

(4-24b) J^{£%^u,,,,sUy) ^ E .9AjdensUy[dsL,dsL])d'x. 

\a\<N-l 

Proof. To prove (|4.23p . we use the deflnition (|4.ip of Un-i and equation p.7cp 
(differentiated with ds) to conclude that 

(4.25) 

l(U^^^^) = 2hcl-l + ^}c.Ul,_,+2e'^'^'{il) EE/ a-")5sA-'^3 
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Inequality (|4.23p now follows from (|4.25p . the definition of Un-i, and the Cauchy- 
Schwarz inequality for integrals. 

(|4.24bj) and (|4.24ap - (|4.24bp follow easily from the definitions (g^Ql), (|4.21a|) - 
(|4.21bp . and the divergence theorem. 

□ 



5. SoBOLEv Estimates 

In this section, we derive Sobolev estimates for the inhomogeneous terms appear- 
ing in equations p.2ap - (j3.2bp . p.3|) . and (|3.7al) - p.7cl) . for the terms in Lemma 
14. 2[ and for the right-hand sides of (|4.23p - (j4.25p . All of the bounds arc derived 
under a smallncss assumption on the norm Sn- These bounds form the backbone 
of Section [71 where they are used to help derive integral inequalities for IAn-i and 
£m, £N;veiocity, and £ N -i;density We coUcct togcthcr all of the estimates in the next 
proposition. 

Remark 5.1. Note that in following proposition, we establish different estimates 
for the special cases = and = 1/3 compared to the other cases. 

Proposition 5.1 (Sobolev Estimates). Assume thatO < < 1/3. Let {L,u^ ,u^) 
be a classical solution to the relativistic Euler equations (|3.2ap - (j3.2b[) on the space- 
time slab [1,T) X R"^, and let Sn be the solution norm defined in Definition \4-l\ 
Assume that S]\[{t) < e on [l,T). Then if e is sufficiently small, the following esti- 
mate for holds on [1, T) : 



e ^ N \velocity , — , 

.^-\n)SN, o<c2<i/3, 

Sn, c2 = 1/3. 

For the time derivatives of the fluid quantities, we have the following estimates 



on 



(5.2b) \\dtL\\H^-.<\ "tc 



(5.2a) \\dtL\\HN-2 <{uj + l)e-^^SN:veloc^ty, cl = 0, 

uj.^-^n)SN + e-"SN, 0<c2<l/3, 
e-"5Ar, cl = 1/3, 

{oj + l)e^^ SN-velocity, = 0, 

(5.2c) |iatu"||H«-i < <( Lo.^-'^in)SN + e~".^-^{n)SN, 0<cl<l/3, 

e-''SN, cl = 1/3, 

{oj -\- l)e^^ SN-velocity, = 0, 

(5.2d) ll^tu^ll^^iv-i < { uje-^^-\n)UN-i + e-'^^SN, < < 1/3, 

c-^Sn, c2 = 1/3. 
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For the inhomogeneous terms u ^^^^a^^ Qabu'^u'' and (3c^ — 2)u)u^u^ from Propo- 
sition [3n[ we have the following estimates on [l,r) in the cases < < 1/3 (the 
case Cg = 1/3 is not directly needed for the remaining estimates): 



CO gabU U 



< 



N, 



(5.3a) 

(5.3b) \\{3cl ~ 2)oju'^u^\\hn < Lje~"SN, 

(5.3c) \\ds{uju"u^) - uju"dsu^\\L2 < ije-^^-\n)SN, {\d\ < N). 

For the inhomogeneous terms A', A"^, and A'-' from Corollarv \3.Sl we have the 
following estimates on [1,T) : 



(5.4a) 
(5.4b) 

(5.4c) 
(5.4d) 



lA'i 



1 < 



uj.^-\n)SN + c^^^Sn, 0<cl< 1/3, 



N, 



1/3, 



{UJ + l)e "^^Sn -.velocity, — 0, 

'-Hn)SN, 0<c2<l/3, 



e-^SN, cl = 1/3, 

(CJ + l)e-'^^'^SN -velocity, 0, 

I A'J'|l^«-i { Lje-^^^-^{n)UN^i + g-^^^Sn, < c2 < 1/3, 



c-'^^^Sn, cl ^ 1/3. 

For the LP' norm of the variation vP = -^Ua'ii"' defined in (j4.8p . we have the 



following estimate on [1,T) : 



(5.5) 



S N -.velocity J2a=l li'^^IU^, 



0, 



^*"IIl^<<^ e"^-i(C!)57vELill^1U^, 0<c^<l/3, 



ci = 1/3. 



For the Lp norms of the inhomogeneous terms ^s., ®% defined in (j4.14b[) . (j4.15p . 
we have the following estimates on [1,T) in the case c^ — : 



(5.6a) 


11^? 






(5.6b) 


i 






(5.6c) 


110 



-2f2N 



<s 



N;velocity I ^ 



~2n 



(0 < |a| < iV- 1), 
(0 < |a| < iV), 



lleslU^ < + l)e-^'''SN-velocity, (0 < < N) . 

Furthermore, we have the following estimates on [1,T) in the cases < < 1/3 
(the case = 1/3 is not directly needed for the remaining estimates): 



(5.7a) 
(5.7b) 



el 

a 

el 



N 



L2 



S 



N 



-20 
-20 



(0 < |«| < iV), 
(0 < \d\ < N). 
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In the case = 0, for the currents J^eiocityl'' '] "'^'^ Jdensityl'^ '] defined in (|4.17ap 
- (|4.17bp . we have the following estimates on [1,T) : 

(5.8) 

E / ^,{jLoc^ty[9du\u',u^),^siu\u',u^)])d'x<iL. + ^)e-■'^^^^^^^ 



\a\<N' 

(5.9) J2 I 5M(^L.s*,[55i,55£])d3x<(a; + l)e-2^^5^. 

In the cases < < 1/3, for the currents J^[-, •] defined in (j4.16p . we have the 
following estimates on [l,r) : 

(5.10) 

Remark 5.2. We have not attempted to derive optimal estimates in Proposition 
15.11 we only derived estimates that are sufficient to prove our main future stability 
theorem. 

Proof. Proofs of (|5.ip - (j5.4dp : We first consider the cases < < 1/3. To prove 
(|5.ip wc use equation (|1.6p . Definition 14.11 Proposition IA-31 (with F{v) ~ a/T+Ij, 
V = gabu'^u'', and v = 0), ProDOsition lA-41 and Sobolev embedding to conclude that 



- < \\9abu''u''\\HN < e^''5ab\\u''\\L.^\\u''\\HN < .^-\n)UN~l < .^-\n)SN, 

where Sab is the Kronecker delta. The estimates (|5.3ap - ()5.3cp and (|5.4ap - (j5.4dp 
follow similarly (use Proposition IA-51 to deduce ()5.3cp ). The estimates ()5.2ap - 
(|5.2dp then follow trivially with the help of equations p.7ap - p.7cp . The cases 
= 0, 1/3 follow similarly. 

Proof of (|5.5p : To prove (|5.5p in the cases < < 1/3, wc use equation (|4.8p . 
(|5.ip . and Sobolev embedding to deduce 



(5.11) 



< 



'^-\n)SN\\u''\\L^. 



The cases = 0, 1/3 follow similarly. 



Proofs of ([5?6a|) - (|5.7bp : We first prove (|5Ja| (in the cases < < 1/3). We 
begin by using equation (|4.14bp to deduce that 



(5.12) 



\4) 



< \\dsh-hs\\^, + \\A^ds[{A'')-'h]-dsh\\^, 



L2 
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See the remarks at the end of Section 12.31 concerning our use of notation for the 
norms of array- valued functions. 

We now deduce the following preliminary estimates for W (L, ^v? , u^)'^ , 

b ^uj(^-^-^gabu''u'', (3c2-2)uOui, (Bc^- 2)^0^2, {3c^^-2)u''u^y , hs uJ{3c^^- 
2)u^(o,^su^,^su'^,^su^^ , |d| < N, A^, and (cxphcit expressions for the 
matrices arc provided in (|3.19ap - p. 20^ *1: 



(5.13) 

l|b|lHiv-i < ujSn 

(5.14) 



\\dsh-hs\\L2 <ojS, 



N 



( ^-\^) \ 



e '-tiff "(r2)5jv 



e^^-\n)SN e^^-^{n)SN e"^-\n)SN\ 



e '"m "(il)5jv 



(5.15) 



/ 1 



^^'^-\n)SN e".^~^in)SN e^l^-\n)SN\ 
N 1 ^-^{n)SN ^-Hn)SN 



1 



N 



(5.16) \mA^)-'\\^.-,<S_ 



1 



(5.17) 
(5.18) 



11^ 

e ^' 



-2n 








N 









N 



\ 






(5.19) 



N 



( e-^ \ 

2f^^-i(f}) e-" 



(SI) e 
\d) e-"/ 
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and analogously for the matrices A^, A^. The above preliminary estimates follow 
from repeated applications of Proposition lA-21 ProDOsition lA-41 the definition (|4.2p 
of Sn, and the estimates ()5.ip - (|5.3cp . 

We now estimate the second term on the right-hand side of (|5.12p . To this end, 
we use Propositions IA-41 and IA-5[ together with the above estimates and Sobolev 
embedding to conclude that for < |q;| < 



(5.20) 



A^ds[{A°)-^h] - dsh < ||A"||l=o * \\diA 



lO\-l I 



*llb| 



N 



\e-^^.ss-^{n)j 

where we write * to indicate that we are performing matrix multiplication on the 
matrices of norms. 

We similarly estimate the third term on the right-hand side of (|5.12p (using 
Proposition Ell] to estimate \\d[{A°)-'^ A^^hn-i < ||(^°)-1||l~ * WOA^'Uhn-i + 
\\d{A°)^^\\HN-i * IIA'^IIl^), thus arriving at the following bound: 

(5.21) 



A'>{{A°y'A-dadsW - ds[{A"r^A-daW] } 

< * {||(A")-l||i=o * WM'Wh--^ + \mA")-^\\HN-i * \\A^L^} * \\daW\\H^ 



N 



-20 



Finally, adding the second estimate in (j5.13p and (|5.20p - (j5.2ip . we deduce 
(|5.7ap . The proofs of (|5.6ap - (j5.6b[) (for = 0) are similar, and we omit the 
details. 

To prove (|5.7b|) . we use equation (|lT5)) . ((5T|) . ([Qc)) - (|5.2dp . and ([E7a|) to 
conclude that 



(5.22) \\&%h-^<gab\\u''d,(^) 



\U II^N + LJ 



(3c2 -2 + 2u°)ua 



l"a|k-!l®5llL^ 



< 



N 



This completes the proof of (|5.7b|) . The proof of ()5.6cp (for ~ 0) is similar, and 
we omit the details. 



Proofs of dSJl) - ((5J0)) : To prove ([5J0)) in the cases < < 1/3, we first recah 
equation ()4.18p . where W = (L, u\u^,u^)^ ds'W = {dsL, d^u^ , dsu^ , dsu^)'^ , 
and as in (|4.8p . =^ -^Uau"' : 
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(5.23) a^(>[w,w]) 



f 2clgab f 9) 



u L + 4c;aj — L 



+ 2{l + ci){ici-l)ujgabu^u'' 



< if < 1/3 



L - 2(1 + + 2(1 + cj)ga6®"it 



Above, the terms 5 = ^5 and © = 0^, are defined in (|4.14bp and (|lJ5)) . We 
now deduce (jS.lOp using the following four steps: i) we ignore the non-positive 
term 2(1 + Cg)(3c^ — l)ujgabU°">J-^ on the right-hand side of (|5.23p : ii) we bound 
each variation L, and 5^5, in L^, and use the estimates (|5.5p . (|5.7ap . and 
(|5.7bp : iii) we bound all of the remaining terms in L°° and use Sobolev embedding 
together with the estimates (|5.ip and (|5.2cp - (|5.2dp : iv) we make repeated use 
of the estimate ||wiW2W3||ii < ||i'i||L°°||'y2||L2||?;3||j;^2, where V2 and V3 are terms 
estimated in step ii), and vi is estimated in step iii). 

The proofs of (|5.8p - (j5.9p (in the case = 0) arc similar but simpler (use the 
expressions (j4.19ap - (j4.19bp ). The proof of (|5.10p in the case = 1/3 is provided 
in the next paragraph. 

Details for the special case ~ 1/3: Wc now address the case = 1/3. Our 
goal is to show that in this case, many of the estimates (|5.ip - (|5.10p are valid 
without the lu ■ ■ ■ terms on the right-hand side. The fact that these estimates 
hold without these terms is a consequence of special algebraic cancellation. This 
cancellation, which is connected to the conformal invariance of the relativistic Eu- 
ler equations when = 1/3, is discussed in Section [9] in more detail. We will 
only demonstrate the cancellation for the estimate (j5.10p : this is the only esti- 
mate needed to derive the crucial inequality (|7.2p below in the case = 1/3. 
In Section ini we will show that when = 1/3, {p,u) verify the relativistic Eu- 
ler equations p.7ap - (jl.7bp corresponding to the metric g if and only if the 

rescaled variables (p', U) (e^^p, e^u) verify the relativistic Euler equations (jl.7ap 

- (jl.7bp corresponding to the Minkowski metric m = e~'^^g. Now if we consider 
the change-of-time-variable ^ = e"^^*-*, T{t = 1) = 1, then (|9.2p shows that we 

have m = -dr^ + Ej=i(f^a;^T- Therefore, it follows that L = In (e4"(*°^V/p), 
jjj e^^'^*°'^'u^ , verify the relativistic Euler equations ()3.2ap - ()3.2bp corresponding 
to the Minkowski metric m = —dr^+J^j^iidx^y and the coordinates {t,x^ ,x^,x^). 
In particular, in these new fluid and time variables, the right-hand sides of (|3.2ap 

- (|5:2b)) vanish, and W = {1 + SabW'U^y/'^ = u^>. Here, Sab is the Kronecker delta 
and the index "r" is meant to emphasize that we are referring to the time com- 
ponent of U relative to the new time coordinate system (while u'^ refers to the 
time component of u relative to the original "i" coordinate system). Furthermore, 
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relative to these new variables, we can define a current analogous to the current 
defined in (|436)) : 

(5.24a) = TTT^^' + 2c'f>'i + (1 + <^lW\-{^^f + ^abU'^uX 

(5.24b) J'^ j^^^L' + 'iclWL + (1 + ci)W[-{U-f + 5abU'U\ 

where tl'^ = e^u^ [/^ = ^SabU^U^ ^ -^gabu^'ii'' = u°. Note that = i°, 
= e^.P, and 

(5.25) dr/^ + da/"" = e'\dtj'' + daJ"-) - e''9^>. 

Setting W {L,U^ ,U^), using the vanishing of the right-hand sides of 
p.2ap - (j3.2bp in the new coordinate + fluid variables, and noting definition (|4.2p . 
our prior proof of (|5.10p (under the assumption that the metric is the standard 
Minkowski metric) yields 

(5.26) V / dr{J^[d3W,dsW]) + da{J''[dsW,dsW])d''x 
Using ((05)) . ((5?26|) . and definition (|42]) . we deduce 



'^x 



(5.27) e'^ V / a^(>[a5W,a5W])d 

= y. I dr{J^[dsyv,dsyv]) + da{J''[doyv,d3yv])d 

3=1 3=1 

which is the desired estimate. 

□ 

6. The Energy Norm vs. Sobolev Norms Comparison Proposition 

In this section, we prove a proposition that compares the coercive properties of 
the norms Un-i, Sn, and Sjq-.veiodty to the coercive properties of the energies fjv, 

^N-jVelocity 7 and S N —l\density • 

Proposition 6.1 (Equivalence of Sobolev Norms and Energy Norms). Let 

N > 3 be an integer, and assume thatSN{t) < £ on [1,T). Then there exist implicit 
constants such that if e is sufficiently small, then the following estimates hold for 
the norms Un~i, Sn, o,nd Sn -velocity defined in Definition \4. 1\ and the energies Ejq, 
£N;veiocity, and £ M -1- density defined in Definition \4.3\ 
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C'2 I c2 ^ c2 „2 _ n 

'^N-velocity N -I- density ~ '-'JV J — 

(72 ^ c2 ^2 = n 

('e 1 '-' N;velocity ^ N ^velocity 

^ + « 52,, 0^< C2 < 1/3, 

Proof. In the cases < < 1/3, to prove Ujf_i + £^ ~ 5^, we only need to show 

3 . 

(6.2) ^ {lISsLlli.+e^^^^lias^i^lli,}^ ^ / j"[dsW,dsW]d^x. 

\S\<N a=l \3\<N-^^^ 

The desired inequahties would then follow easily from (|6.2p and Definitions 14.11 
and 14.31 To prove (|6.2p . for notational convenience, we set W (L, u^, m^, u'^)-^, 
■uJ =^ SsM^, L =' and as in ()4.8p . we define ■^Uau"'- We now recall 

definition 

(6.3) j° = j^^^L^ + 2c2itOL + (1 + c2)^05„^u"it'^. 

Using < < -i, (|5.ip . Sobolev embedding, and the simple estimate |2Cgu"i| < 
2{ilci) i^^ + 2cs(l + c2)(^J^MaM°) , it follows that 

(6.4) j" > ^M_1/^L' + (1 + c^)^..^'^^'' - (1 + cl){u" + 2c^)(l.„u")' 

A reverse inequality can similarly be shown. Integrating these inequalities over 
M'^ and summing over all derivatives d^W with |cJ| < A*", we deduce (|6.2p . The 
remaining estimates in (|6.ip (in the cases = 0, 1/3) can be proved similarly. 

□ 

7. Integral Inequalities for the Energy and Norms 

In this section, we use the Sobolev estimates of Section [5] to derive energy and 
norm integral inequalities. These inequalities form the backbone of our future 
stability proof. 

Proposition 7.1 (Integral Inequalities). Let {L,u^ ,u^) be a classical solu- 
tion to the relativistic Euler equations (j3.2ap - (|3.2bp on the spacetime slab [1,T) x 
M.^. Let Un-1-, Sn, and SN;veiocity be the norms in Definition \4-l\ and let En, 
£N;veiocity, o,nd £ N -i;density bc the energies in Definition \4-3\ Let N > 3 be an 
integer, and assume that 5Ar(t) < e on [1,T). //c^ = and e is sufficiently small, 
then the following integral inequalities are verified for \ < ti < t < T : 

(7.1a) El,,,i,,,ty{t) < £'N.,,eiocUyiti) + C f {oj{s) + 1 } e-^^f^)^^,.,,,^,,,^ (s) ds, 



(7.1b) £^_i^,,„,,,,(t) < £],^,.,dens^ty{tl) + C / {uj{s) + l} e-'''^^^ S^^) ds . 

J S — ti 
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// < < 1/3 and e is sufficiently small, then the following integral inequality 
is verified for 1 < ti < t < T : 

(7.2) £%{t) < £%{t,) + C [' w{s)£S-\n{s))Sl{s)+er''^'^Sl{s)ds. 

absent if — 1/3 

//O < < 1/3 and e is sufficiently small, then the following integral inequality 
is verified for 1 < ti < t < T : 

< for large 



(7.3) Ul_,{t)<Ul_,{t,)+2f^_^ {3c2-l + fflM},(,)Z^^_^(,) 



ds 



C ^* 



• / u{s).^^^{n{s))Ul_^{s) + e-^^^"'^ ^[n[s))S%[s) ds. 

J s=ti 

Proof. To prove ([7^ . we first use (|T23| and (|5.4dp to deduce that 



(7.4) 



< for large 



^^{Ul^,) < 2 {3cl - 1 + :|M }^u'^_^ + 2e^m)U^^, || A-||^. 
< for large 



Integrating (|7.4p from ti to t, wc deduce (|7.3p . 

Inequalities ((7?Ta)) - (j7.1b|) and follow similarly from (|4.24ap - (|4.24bp . 
and jnil) - dSH), ((5T01) . 

□ 

8. The Future Stability Theorem 

In this section, we state and prove our future stability theorem. The proof is 
through a standard bootstrap argument based on the continuation principle (Propo- 
sition [3]4]) and the integral inequalities of Proposition [TT] 



Theorem 8.1 (Future Stability of the Exphcit Fluid Solutions). Assume 
that < < 1/3, that the metric scale factor verifies the hu'Dotheses \Al\ - \A3\ from 

Section[l\ and that N > S is an integer. Let L = L\t=i == hi (^p / pj \t=i 
(j — 1,2,3) be initial data for the relativistic Euler equations (j3.2ap - p.2b[) . and 
let Sf^it) be the norm defined in Definition \4-l\ There exist a small constant eg 
with < eo < 1 and a large constant C* such that if e < eg andSjsii^) — C^^e, then 
the classical solution (L = ln(e3(i+'==)f2p/p),u) provided by Theorem \ 3. 1\ exists for 
{t^x^ ,x'^ ,x^) e [l,oo)xM'^ and furthermore, the following estimate holds fort > 1 : 



(8.1) 



Swit) < e. 
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In addition, Tmax = oo, where Tmax is the time from the hypotheses of Proposition 

[M 

Proof. Wc provide full details in the cases < < 1/3; the remaining case ~ 
is similar and simpler, thanks to the fact the evolution of the decouples from 
that of L. By Theorem 13. 1[ if C* is sufficiently large and e is sufficiently small, 
then there exists a non-trivial spacetime slab [1, Tiocai) x upon which a classical 
solution exists and upon which the following estimate holds: 

(8.2) SNit) < e. 

We define 

(8.3) 

T =* sup {Tiocai > 1 I The solution exists on [l,Tiocai) and (|8?2)) holds on [1, Tiocai)}- 

We will show that if C* is sufficiently large and e is sufficiently small, then T ~ oo. 
Throughout the remainder of the proof, we use the notation 

(8.4) e = 5w(l) = C-ie. 

We first use the bootstrap assumption (j8.2p (under the assumption that e is 
sufficiently small), Proposition l6.H the inequalities (j7.2|) - (|7.3p . and the hypotheses 
lAll - IA3I from Section [T] to derive the following inequality, which is valid for all 
sufficiently large ti and t verifying ti < t < T : 

(8.5) Sjf{t)<CSNiti) + Ce^ f uj£§'\n) +e-^l!^{n) ds. 

Js = ti V ' 

absent if = 1/3 

Now our hypotheses lAll - IA3I on imply that the integrand on the right-hand 
side of (|8.5p is integrable in s over the interval s G [l,oo). Therefore, if ii < T and 
ti is sufficiently large, it follows that 

(8.6) Slit) <CSNiti) + ^e^, te[t,,T). 

It remains to derive a suitable bound for 5Ar(ii); this will be a standard local- 
existence- type estimate. To this end, for this fixed value of ii, we again use the 
bootstrap assumption (|8.2p . Proposition 16. 11 and (|7.2p - (|7.3p to derive 

(8.7) Sjfit) < + c{ti) [ S%{s)ds, te[l,t,). 

Js=l 

Applying Gronwall's inequality to (|8.7p . we deduce that 

(8.8) S%{t) < Cfe'^''^\ t€[l,t,). 

Furthermore, we note that by Proposition 13.41 Sobolev embedding, and the conti- 
nuity of SNit), it follows from (|8.8p and definition (|8.4p that if e is sufficiently small 
and C, is sufficiently large, then 
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(8.9) <i < T. 

Combining (|8.6p and (jS.Sp and using definition (|8.4p . we deduce that 

(8.10) 5^(t)<Ce^(*i)^+l,2^ ie[l,T). 

Therefore, if e is sufficiently smaU and is sufficiently large, it follows from (|8.10p 
that 



(8.11) Sl{t)<e^, te[l,T). 

Note that (|8.1ip is a strict improvement over the bootstrap assumption (|8.2I) . Thus, 
using Proposition 13. 4[ Sobolev embedding, and the continuity of SN{t), it follows 
that r = cx) and that ([5TT|) holds for i e [1, oo). 

The case = can be handled similarly using inequalities (|7.1ap - (|7.1b[) . 

□ 

9. The Sharpness of the Hypotheses for the Radiation Equation of 

State 

In this section, we show that if = 1/3, hypothesis [Al] from Section [T] holds, 
and e"^^''^ ds = oo, then the background solution p = pe~'^^, = 6q to the 
relativistic Euler equations (|1.7ap - (jl.Tbp on the spacetime- with-boundary ([1, oo) x 
M'^,g) is nonlinearly unstable. The main idea of the proof is to use the special 
conformally invariant structure of the fluid equations when = 1/3 to reduce the 
problem to the case in which the spacetime is Minkowskian; we can then quote 
Christodoulou's results |Chr07j to deduce the instability. We begin by noting a 
standard result: that the change of time variable 

(9.1) ^ = e-^^(*) 

^ ' dt 

allows us to write the metric ()1.2p in the following form: 



(9.2) g = e2"(*°-) ( - dr^ + ^(dxJ)^) = e2^^(*°-)m, 

where m — —dr^ + X]^=j('^^"')^ the Minkowski metric on [l,oo) x equipped 
with standard rectangular coordinates. The following elementary but important 
observation play a fundamental role in the discussion in this section: if the hypoth- 
esis lAll is verified but J"^^ e"^^*^"' ds = oo, then after translating r by a constant if 
necessary, it follows that 



(9.3) r : [l,oo) [l,oo), t ^ T{t) 

is an autodiffeomorphism of [1, oo). 

We now prove the conformal invariance of the fluid equations when = 1/3; 
this is a standard result, and we provide the short proof only for the convenience 
of the reader. 
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Proposition 9.1 (Conformal Invariance of the Relativistic Euler Equa- 
tions When = 1/3). Let g he the conformally flat metric defined in (j9.2p . Then 
{p, u) is a solution to the relativistic Euler equations p.7ap - (jl.7bp corresponding 
to the metric g if and only if (p', U) is a solution to the relativistic Euler equations 
()1.7ap - (|1.7bp corresponding to the Minkowski metric m. Here, the rescaled fluid 
variables (p',U) are defined by 



(9.4) P = e p, U ~ e u. 

Proof. Throughout this proof, V denotes the Levi-Civita connection corresponding 
to the Minlcowslci metric m and D denotes the Levi-Civita connection corresponding 
to g = e^^m. As discussed in the beginning of the article, the relativistic Euler 
equations (|1.7aP - (|1.7bp are equivalent to the equations DaT"'^ = plus the 
normalization condition gapu°'u^ ~ —1 (see (|1.3p and (jl.Sp V In the case p = 
l/3c^p, we note that 

(9.5) T^'' = ^pu'^u'' + \p{g-'r - Iru^W" + ipe-2^(m-l)^^ 

and we define the following rescaled "Minkowskian" energy-momentum tensor: 

(9.6) Tl^^^ = e^'^T^"' = ^p'VU'' + ip'(m-i)^'', 
where p' and are defined in (|9.4p . We also note that by (|9.2p . 

(9.7) mc.pU'^U^ = -1. 

The key step is the following identity, whose simple proof we omit: 



(9.8) DaT"'' = e-*^^^Va(e^*^^T"^) = er^^V aT^J'^y 

We remark that the proof of the identity (j9.8p heavily leans on the fact that 
gafsT"^ = for the equation of state p ~ (l/3)p. It now follows from (|9.8p that 



(9.9) D^T"'' ^ ^^t;^J^) = 0- 

According to the remarks made at the beginning of the proof, this completes the 
proof of the proposition. 

□ 

In the next theorem, we will recall some important aspects of Christodoulou's 
shock formation result jChr07j . To this end, we need to introduce some notation. 
In order to emphasize the connections between Proposition 19. 11 Theorem 19.11 and 
Corollarv 19.21 we will denote the energy density by p', the four velocity by U, and 
the spacetime coordinates by (r, 3). Let (p',^7^^7^^73) be initial data (at 

time T = 1) for the relativistic Euler equations (|1.7ap - (jl.7bp on the spacetimc- 
with-boundary ([l,oo) x M^,m), where m = — dr^ + X]j=i('^^"')^ standard 
Minkowski metric on M^. Let p > be a constant background density. Let Br C 
denote a solid ball of radius | < r < 1 centered at the origin in the Cauchy 
hypersurface {(r, x^, x^, x'^) | r = 1} ~ K'^ (embedded in Minkowski spacetime) 
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and let dBr denote its boundary. We define the order-1 Sobolev norm SBi\Br ^ 
of the data over the annular region in between _B,. and the unit ball as follows: 

3 

(9.10) Sb,\b,. +'^\\U^\hhb,\b,.)- 
We define the standard Sobolev norm of the data as follows: 

3 

(9.11) Vm = \\p'-p\\h^'+Y.\\^'\\h^'- 

We also define the following surface + annular integrals of the data: 
(9.12) 

Q{r) [ r\{p' -p) + -^pU'nA da + [ 2{p' ~ p) + -^pWN, d'x. 

Above, r =' denotes the standard radial coordinate on ffi.^, and N 

denotes the outward unit normal to dBr- 

Theorem 9.1 f |Chr07j Theorem 14.2, pg. 925). Consider the relativistic Euler 
equations (|1.7ap - (|1.7bp on the manifold-with-boundary [l,oo) x equipped with 
the Minkowski metric and rectangular coordinates (t, ) : m = —dr + 

X]j=i('^^"')^- p > be a constant background density. Then the explicit solution 

p' = p,{U'^ ,1}"^ ,U^) = (1,0,0,0) is unstable. More specifically, there exists 
an open family of arbitrarily small ( non-zero ) smooth perturbations of the explicit 
solution's data which launch solutions that form shocks in finite time. Here, 
denotes the time component of U relative to the coordinate system {t,x^ ,x'^ ,x'^). 

Even more specifically, there exists a large integer M, a large constant C > 0, 
and a small constant e > such that the following three conditions on the initial 
data together guarantee finite-time shock formation in the corresponding solution: 



(9.13a) Vm < e, 

(9.13b) Q{r)>C^/v^{^V^+VT^}S(^B^\B^), 
2 

(9.13c) ^<r <1. 



Furthermore, there exists a uniform constant C > such that the above three 
conditions guarantee that a shock forms in the solution before the time 

(9.14) Tmax{r) = exp 

Remark 9.1. Note that if we are given any data such that Q{r) > 0, then if 
we rescale its amplitude (more precisely, the amplitude of its deviation from the 
background constant solution) by a sufficiently small constant factor, then (|9.13ap 
and (j9.13bp are both verified by the rescaled data. This follows from the fact that 



C'il-r) \ 

\ Q{r) r 
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for a fixed r, I?j\f , 5(Si\_b^), and Q{r) shrink linearly with the scaling factor, while 
the right-hand side of (j9.13bp shrinks like the scaling factor to the power 3/2. 

The following corollary follows easily from (|9.1|) - (|9.2p , Proposition 19.11 and 
Theorem [9Tl 

Corollary 9.2 (Nonlinear Instability of the Explicit Solutions When = 1/3). 

Assume that = 1/3. Consider the relativistic Euler equations (|1.7a|) - ()1.7bp on 
the manifold-with-boundary [l,oo) x MP equipped with a Lorentzian metric g = 
~dt^ + e2f2(*) Ej=i(t^a;')', f^(l) - (as m (|1.2|) ). Assume that hypothesis [Al\ of 
Section[l\ holds, and that J^i^^^^'^^ ds = oo. Then the explicit solution (|1.8p is 
unstable. More specifically, there exist arbitrarily small (non-zero) perturbations of 
the explicit solution's initial data which launch perturbed solutions that form shocks 
in finite time. More precisely, if the open conditions (|9.13a[) - (j9.13cp are verified 
by the data, then a shock will form before the conformal t— coordinate time 

(9.15) Tr,iax{r) = exp 

In the original time coordinate t, the shock will form before tmaxif), which is im- 
plicitly determined in terms of Tmaxi^) via the relation 

(9.16) w(r) 

Remark 9.2. Unlike the proof of Theorem 18. 1[ the proof of Corollary 19.21 is not 
easily seen to be stable under small perturbations of the metric g. In order to show 
that the proof is stable, one would need to show that Christodoulou's proof of 
Theorem 19. H is stable under small perturbations of the Minkowski metric. 
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Appendices 

A. SOBOLEV-MOSER INEQUALITIES 

In this Appendix, we provide some Sobolev-Moser estimates. The propositions 
and corollaries stated below are standard results that can be proved using meth- 
ods similar to those used in |H6r97[ Chapter 6] and in the Appendix of [KMSlj . 
Throughout this Appendix, LP = LP{R^) and ff*^ = H'^^R^). 

Proposition A-1. Let AI > be an integer. If {va}i<a<i oltc functions such that 
Va G L°° ,\\d}'^^\a\\L'^ < OO for 1 < a < I, and (5?i,--- ,(5; are spatial derivative 
multi-indices with \ai \ + ■ ■ ■ + \ai\ = M, then 



C'{l~r) \ 

\ Q{r) r 
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(A.l) \\{da,v,){ds,v2) ■ ■ ■ {da,vi)\\L^. < C{l,M) ^ (\\d^''\k.\\L- J] I^^IU-) • 

a— 1 h^a 

Proposition A-2. Let M > I be an integer, let M. be a compact set, and let 
F e C^^(^) be a Junction. Assume that v is a function such that viM?) C ^ and 
dv e H'^'-\ Then d{F o w) e H^^-i, and 

M 

(A.2) ||9(Fo^;)||^„-, <C(M)||a^;||HA.-i5^|i^(')Ulkll^~i. 

1=1 

Proposition A-3. Let M > 1 be an integer, let A be a compact, convex set, and 
let F G C^^ (A) be a function. Assume that v is a function such that u(R'^) C M. 
and V — V E H'^^ , where v E ^ is a constant. Then Fov — FovE H^^ , and 

(A.3) 

M 

\\Fov -FovWhm < C(A/){|f(i)U|K> - v\\l2 + \\dv\\m'-^ E 

1=1 

Proposition A-4. Let M > 1,1 > 2 be integers. Suppose that {va}i<a<i are 
functions such that Va G L°° for 1 < a < I, that vi E H'^^ , and that dva £ H^^^^ 
for 1 < a < I — 1. Then 

(A.4) 

i-i i-i 

a=l a=l b^a 

Remark A.l. The significance of this proposition is that only one of the functions, 
namely vi, is estimated in L^. 

Proposition A-5. Let M > 1 be an integer, let ^ be a compact, convex set, and let 
F G C^^(^) be a function. Assume that vi is a function such that vi{M.^) C A, that 
dvi e L°°, and that d^^^'>vi € Assume that V2 € L°° , that d^^''^\2 e L^, and 
let d be a spatial derivative multi-index with with \d\ — M. Then ds {{F o Vi)v2) — 
{F o vi)dsV2 G L^, and 

\\dsi{Fov,)v2)-iFovi)dsV2\\L- 
(A.5) 

M 

< C{M){\F^'^U\\dvi\\L^m^''~'^V2\\L^ + \\v2\\L^\\dv^\\HM-^ ^ \F^''> U\\vi\\'l2} . 

1=1 
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